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Abstract 

The effects of quantum and thermal corrections on the dynamics of a damped nonlinearly kicked 
harmonic oscillator are studied. This is done via the Quantum Langevin Equation formalism 
working on a truncated moment expansion of the density matrix of the system. We find that the 
type of bifurcations present in the system change upon quantization and that chaotic behavior 
appears for values of the nonlinear parameter that are far below the chaotic threshold for the 
classical model. Upon increase of temperature or Planck's constant, bifurcation points and chaotic 
thresholds are shifted towards lower values of the nonlinear parameter. There is also an anomalous 
reverse behavior for low values of the cutoff frequency. 

PACS numbers: 05.45. +b, 03.65. Sq, 05.40. +j, 42.50.Lc 
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I. INTRODUCTION 



With the impressive results obtained in the theory of classical chaos, the interest in the 
manifestations of chaos in quantum systems has grown over the past years. Although there 
is not a clear difference between chaotic and regular dynamics for quantized Hamiltonian 
dynamics, some insight has been gained with the quantization of conservative models with 
chaotic classical counterparts. Among the peculiar phenomenology that has been found we 
have the appearance of scars, which are wavefunctions with enhanced amplitude along the 
paths of classical unstable periodic orbits l|, SSQ]- Further connection between classically 
unstable orbits and and quantum behavior has been given by Gutzwiller's trace formula for 
the den.it^ ot .tate.]! which ha. been appUed with some success to many l^yperbolic 
systems jTi]. A connection has also been found between the nature of the classical dynamics, 
be it regular or chaotic, and the statistics of the energy spectrum, which obeys the general 



distributions found in Random Matrix Theory 

However, for systems that include dissipation, not much work has been done. This 
is an important subject of study because, besides its theoretical interest, it is of central 
importance in several interesting applications (SQUIDS, tunneling phenomena, etc.) The 
main problem with such systems is the presence of the dissipative term, which makes the 
quantization process not feasible by standard methods. Approaches to this problem have 
been proposed by many authors over the last three decades. Kostin proposes the use of a 
nonlinear Schrodinger equation for the process of quantization , but it has some problems 
with the superposition principle and yields results of an unphysical nature. Dekker ^| uses 
a canonical quantization method based on complex variables, with some degree of success, 
but the physical basis of his theory is not very clear. Maximal Entropy Methods 1^ have 



also been proposed, using a coupling between quantum and classical degrees of freedom. 
This approach ignores the effects of thermally driven fluctuations. 

In general, the only fully successful and physically sound quantization scheme for dissi- 
pative systems is the so-called system-plus-reservoir approach, in which one considers the 
system coupled to a large reservoir or heath bath. Since the entire system plus reservoir 
may be considered as closed, one can apply to it the standard rules of quantization, paying 
of course the penalty of dealing with an infinite number of degrees of freedom. Dissipation 
comes about by the net transfer of energy from the system to the reservoir, and ther- 
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mal fluctuations appear as a result of the random energy exchange from/to the heat bath. 
These approaches ensure the non-violation of commutator rules, and include the appropriate 
fluct uat ion-dissipat ion relations . 

Some common approaches using the system-plus-reservoir picture are the influence- 
functional methods introduced by Feynman and Vernon Il3| , that have been used by Caldeira 

n n 

and Legget to study quantum Brownian motion [14| and dissipative tunneling effects jl5| : 
generalized master equa.tiqns [l^, used to describe damping phenomena in quantum optics 



and spin dynamics 17, [l^; and Quantum Langevin Equations used by Ford et al. to 
study dissipative quantum tunneling at zero temperature 

In the area of nonlinear (and classically chaotic) dissipative models, there are some results 
that point to the conclusion that quantum effects act as sources of noise. As a result, 
it is expected that the threshold for chaos be lowered in a continuous fashion upon the 



2l| , who worked 



introduction of quantum corrections. This was reported by Graham and Tel 
in the master equation formalism to quantize Henon's map, and by Cerdeira et al. 

[2^1, who built a quantum mapping of Wigner functions for a kicked oscillator. More recently, 
Cohen [2^ studied the quantum dissipative kicked rotator using influence functionals, while 
Dittrich and Graham used the master equation formalism for the same purpose j25|. Liu 



and Schieve 



26l | use the Quantum Langevin Equation approach to show that a chaotic 



attractor is observed in a quantum system whose classical counterpart is non-chaotic. In 
contradiction to these works, there were some results showing an increase of order (a raising 

n 

of the chaoticity thresholds) for a kicked dissipative oscillator |27l |. 

In this work we study the behavior of a periodically kicked harmonic oscillator, with 
dissipation, and in the presence of quantum and thermal fluctuations. Our interest is mostly 
to find how the bifurcation points and chaotic thresholds in parameter space change when 
increasingly large thermal and quantum effects are included. To do this, we generate a 
mapping of expectation values, using a quantum Langevin approach, that has an immediate 
analog with its classical limit. This Expectation Values Mapping (EVM) constitutes in 
principle a complete description of the quantum dynamics, based in a discrete set of variables, 
and has the advantage of exhibiting a clear hierarchy in its terms, introducing from the 
beginning mean values and a sequence of moments of the fluctuations (the completeness of 
the description given by the EVM will be shown in appendix A). 

In section II we introduce the classical model, a harmonic oscillator with a periodic "kick" 
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from a nonlinear potential. In section III we will use the Langevin equation formalism 
to obtain a quantal version of the classical map of the preceding section. The resulting 
(truncated) EVM will include both thermal and quantum corrections; the change in the 
behavior of the dynamics as function of ksT and h is examined. In section IV and V we 
show and discuss the results obtained and give some concluding remarks. 

II. CLASSICAL KICKED HARMONIC OSCILLATOR WITH DAMPING 

Starting with the Hamiltonian for a kicked oscillator 

H ^ ^ + -mn''x^ + V{x)y^S{t-nT), (2.1) 
2m 2 ^ 

where r is the period of the kick, we introduce the dissipative force 270; to obtain the 
equations of motion in the form 



p = -mQ^x - 27^ - V'{x) ^ 6{t - tit), (2.2) 



X ^ ^. (2.3) 
m 

Prom now on we will set the units of mass and time so that m — 1 and D, — 1. For times 
t ^ riT, the kick has no effect so we get the equation 

x + 2-fx + x^O, (2.4) 

whose well known underdamped (7 < 1) solution is given by 

x{t) = e-'^\Acosu;ot + Bsinuot), (2.5) 

p{t) — e"'"'* [u;o{—A sin coot + B cos ujot) — ^{A cos coot + B sin LOot)] , (2.6) 

where cuq — -\/r^^7^, and A and B are given in terms of initial conditions. Evaluating these 
constants, and defining an evolution matrix T by 

Tqq = exp(-7T)(cos(a;oT) +7sin(a;oT)/u;o), (2.7) 

Tqp = exp(-7r) sin(a;or)/a;o, (2.8) 

TpQ = -exp(-7r)(u;o + 7^/'^o)sin(a;or), (2.9) 

Tpp = exp(-7T)(cos((x;oT) - 7sin(c^oT)/c^o), (2.10) 
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we can write a discrete linear map corresponding to Eqs. ()2.51 

^n+l = TqqX^+TqppI, (2.11) 

Vn+1 = TpQX^ + Tppp^, (2.12) 

where we have set as usual x~ = x{nT — O"*"), x^ = x{nT + O"*"), with similar definitions for 
p. We need now to calculate how the kick affects the dynamics, in order to complete the 
map. This can be done by integrating Eq. ()2.2|) from just before to just after t = {n + l)r. 
Over this region x is continuous and thus we obtain 

Pt+l-Pn+l = -y'i^n+l), (2-13) 

Successive application of the mappings ()2.5^ 12. (ij) and ()2.1H^ I2.14j) will give us a map that 
yields the values for x and p just after each kick (dropping from now on the superscripts): 

Xn+l = TqQXn+TqpPn, (2.15) 
Pn+l = TpQ X„, + Tpppn - V'{Xn+l). (2.16) 

We have chosen the force acting on the kicks as 

V'{x) = -K)e-("/'')', (2.17) 

which gives a potential bounded from below. From now on we will set the units of length so 
that a = 1. This force is maximum at x = 0, so that the closer the harmonic oscillator gets 
to its rest position, the stronger the repulsive effect of the kick becomes. We are interested in 
a region of parameter space such that 7 -C tu, and where clear chaotic behavior can be found. 
After some search, we have chosen 7 = 0.03 and r = 10.0. For bifurcation diagrams given 
in terms of Vq, we find a well defined transition to chaos around Vq ~ 5.69287. It should be 
noticed that there is some coexistence of attractors in this sector of parameter space; these 
attractors evolve from order to chaos following Feigenbaum scenario (see Fig. (1)). For the 
fully chaotic region, the attractors are Henon-like (see Fig. (2)). 

The later map will be quantized in the next section in order to see the temperature and 
quantum effects on the dynamics. 
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III. QUANTIZATION OF THE KICKED HARMONIC OSCILLATOR WITH 
DAMPING 



We use the Langevin equation formalism for the process of quantization. This equation 
is obtained from the quantization in the Heisenberg picture of a closed system given by a 
relevant degree of freedom Q, in contact with a large (in the limit, infinite) bath of harmonic 
oscillators Qj. It is given by 

Q + 2 r r(t - t')Qdt' + W'{Q) = F{t), (3.1) 

where Q{t) and F{t) are operators, W{Q) is the potential and F{t) is a random force that, 
together with the memory function T{t), depends on the characteristics of the bath. These 
oscillators are assumed to be in thermal equilibrium at some initial time tg- The memory 
'unction T{t), and the expectation value and symmetric correlations of F{t) are given by 
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{Fit)) = 0, (3.2) 
{F{t)F{t')+ F{t')F{t)) = E^^.^fcoth(^)cosh(t-t')], (3.3) 

r{t) = ^EmjCu|cos(t<;jt)e(t). (3.4) 

where rrij and Uj are the mass and frequency of the j'th oscillator. These expressions imply 
a quantum fluctuation-dissipation theorem. It is assumed that no correlation exists between 
F and Q or P. This implies ignoring some switch-on processes, which arise due to our 
assumption that the bath per se is in thermal equilibrium at time tg- 

From now on, we will use a bath of equal masses and a white {i. e., tu-independent), 
and continuous frequency distribution, which is the case for constant friction or Markovian 
approximation. In this case, 

T{t) = 2j6{t), (3.5) 

{F{t)F{t') + F{t')F{t)) = ^ r dujhujcoth { 4^ ) cos[cu(t-t')], (3.6) 

vr Jo ylKBi j 

and Eq. 1)3.11) reduces to the simpler form 

Q + 2^Q + W'{Q) = F{t). (3.7) 
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The drawback of using the Markovian approximation is that it will generate some divergences 
since we now have contributions from arbitrarily high frequencies in the bath. We will use 
the same harmonic oscillator-plus-kick potential indicated in the Hamiltonian (j2.1|) and 
Eqn. (j2.17|) . Between kicks, the (underdamped) solution of Eq. (j3.7|) is given by 

Q{t) = e-^*(Acoswot + -BsincJot) + /(t,cuo,7), (3-8) 
P(t) = Q(t), (3.9) 

where the solution for the non-homogeneous part has been obtained from a Riccatti equation 
and we have defined ujq in an analogous way as for the classical case. The operators A and B 
are given in terms of initial conditions, evaluated at some initial time to (not to be confused 
with te), and the inhomogeneous term is 

/(t,cJo,7) = — [' dt'e-^^'-''hmiUo{t ~ t')F{t'). (3.10) 

UJQ J to 

For the moment P we get 



P{t) = -76"^* [Acos(cuot) + -Bsin(cuot)] 



+uJoe 



-7t 



A sin(cjot) + B cos(co'ot)] + g{t, cuq, 7), 



(3.11) 
(3.12) 



where 



1 /■* 

g{t,uJo,l) = — / dt'e~^''*~*'^ [--fsmuJo{t-t') +uJocosuJo{t-t')]F{t'). (3.13) 

UJq J to 

Now we will take the integration limits to be to ~^ and t ^ {n + 1)t~, so we can 
evaluate the initial conditions and get 



Qn+l 

p: 



Tqq + Tqp Pt + f{T, 000,1), (3.14) 

„+i - TpQQ+ + TppP+ + g{T,u;o,i), (3.15) 

with the same matrix elements Tij defined before. Defining the fluctuations as 

6X = X - {X), (3.16) 



we can now write down the evolution equations for expectation values and fluctuations: 

[Tpq Tpp) [P+ ) 
Tqq Tqp ] I 6Ql \ I f{T,uJo,i) 




TpQ Tpp 



5Pt 



(3.18) 



since both / and g have zero expectation value. Hereafter we will use non-bold letters to 
denote expectation values. In order to find what happens over the kick, we note that Q is 
continuous and P has an increment, so 

Qn+l ~ Qn+1 ~ Qn+li (3.19) 

PUi = P-+i-V\Q,,^,). (3.20) 

Since the potential term has an operator argument we will perform a second order Taylor 
expansion. This is the first truncation we need to carry out in order to close the system of 
equations. From this we get 

V"'(n 1 

PUl = Pn-,1 - V\Qrr+^) - V^"(Q„+l)5Q„+i ^^(^Qn+l)'- (3-21) 

The potential and its derivatives can now be evaluated numerically since we have no longer 
an operator as argument. However in order to have an EVM suitable for numerical work, as 
in Eq. ()3.17j) . we need to take expectation values of all operators. In order to do this we have 
to calculate, among others, the map for the fiuctuation {SQn+iY so that we can evaluate the 
above equation for Pn+i- We proceed by evaluating (^Qn+i)^ and made another truncation 
to second order in 6, so that we get a map which is also second order. In an analogous fashion, 
we obtain the mappings for the fiuctuations (5P^i)^ and (5P^i 6Qn+i + SQn+i SPn+i) and 
after some rearrangement of terms we have the 5 expectation value maps for Q and P and 
for the fiuctuations (some details are given in Appendix B) 



Qn+i — TqqQn + Tqp P^, (3.22) 
= TpQ g„ + Tpp P+ - V'iQr,^,) - ^'"^^^"^'^ (^g.+iOg.23) 

{SQn+ir = T^Q {5Qnf + T^p i5P;ir (3.24) 

+Tqq Tqp {5Qn 5P^ + SP^^ 5Qn) + f-f, 

{SP;l^,f = RIq i5Qrf + Rip (5P+)^ (3.25) 
+RPQ Rpp {6Qn 6P^ + 5P+ 5Qn) + K+i 

6P^_,, 5Qn+i + 5Q„+i = 2RpQ Tqq (5g„)' + 2Tqp Rpp {5P+f (3.26) 

HRpqTqp + RppTQQ){5Qn6P:[ + 5P+5Q„) 

il'n+l)- 

Here the following definitions have been made: 

hn-f + f-K = {K-f + f-K), (3.27) 
8 



h 



■n 



9 - V"{Qn)f 



(3.28) 



RpQ — TpQ — V"{Qn+l)TQQ, 

Rpp = Tpp — V"{Qn+i)TQp, 



(3.29) 



(3.30) 



and cubic and higher order combinations of /i„, /, 6P and 6Q have been neglected. In 
principle, a complete description of the quantum map could be given if one could use all 
moments {PmQm)- As an extra approximation, a cutoff frequency has to be introduced in 
order to calculate expressions which depend on Eq. ()3.3p . where the divergences introduced 
by the use of the Markovian approximation show up. We will now use this set of equations to 
investigate the effects of quantum and thermal contributions to the dynamics of the system. 
The parameters and potential used are the same of the classical map. 

IV. RESULTS 

In this section we present the results obtained using the maps developed in the previous 
section. We analyze the shifts that occur in the threshold of chaos for the system as function 
of the temperature ksT and h. We also investigate how the cutoff frequency influences this 
behavior. 

In Fig. (3) a bifurcation diagram for small values of h and fc^T is shown (small compared 
with the action and energy scales of the problem). There is, as before, some coexistence 
of attractors. These are included using simulation of the map over many different initial 
conditions. As expected, large segments of this diagram are almost equal to the classical 
counterpart (Fig. (1)). There are however, four very important differences: first of all, 
and quite visible in the figure, there are chaotic regions that appear far below the classical 
threshold of chaos (Vq = 5.69287). One of them, at Vq ~ 4, is associated to the second 
bifurcation of the classical diagram. It shows a loss of stability for the quantum system in a 
parameter sector close to a point where the classical system has a zero Lyapunov exponent. 
This behavior is shown in Figs. (4), which shows the bifurcation diagram, and (5), which 
compares the largest Lyapunov exponent for the quantum and classical maps in this region. 
Another appears for Vq around 2.5 and seems to terminate a stable fixed point branch. A 
similar one appears for Vq ~ 5.55). 

Second, and also apparent in Figs. (4) and (5), is the appearance of a Hopf bifurcation 
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just before the chaotic threshold for the quantum map. These transitions are not allowed 
for the classical system, since is dissipative and two-dimensional. Here however, they are 
preferred over the standard period doubling cascade. We have verified the appearance of 
simple non-chaotic cycle attractors, typical of Hopf bifurcations, for this transition. It should 
be noticed that Hopf bifurcations are also present for all other transitions from periodicity 
to chaos in this map, of which four are visible in Fig. (3) ( Vq ^ 2.5, Vq 4, fti 5.55 and 
Vo ^ 5.7.) 

The third main difference between the two bifurcation diagrams is given by the fact that, 
although the system is still dissipative, there seems to be no clear bound for its variables. 

Here we are showing Q only for a finite interval, but the numerical simulations generate 
values of Q as high as several thousands. This peculiarity affects some chaotic regions but 
not all of them. 

Finally, there is no evidence of periodic windows or band mergings inside the chaotic 
region. There seems to be no topological changes in the chaotic attractor once the last 
threshold is crossed. Also, a great loss of detail can be clearly inferred from the bifurcation 
diagram. This is of course even more evident in the attractors themselves, as seen in the 
example given in Fig. (6). This attractor (at least its projection on the Q — P plane) still 
resembles the typical attractor for the classical map. Here we are also cropping the figure. 

Most of these differences are probably associated with the increase in dimensionality of 
the map, resulting from the use of the EVM. As mentioned, a Hopf bifurcation would not 
be possible in the classical system. The behavior of the quantum map, for higher values of 
h and/or ksT, is similar to the one already described, except for an expected growth of the 
chaotic strips. 

It is also worth mentioning that the attractor (at least its projection on the Q — P plane) 
still resembles the typical attractor for the classical map (see Fig. (6)). It of course becomes 
diffused in the plane, since it is just a projection. Notice that we are strongly cropping this 
figure, since values of P and Q spread up to the thousandths. 

The changes we have mentioned up to now signal qualitative differences between classical 
and quantum maps. Besides these, we have also studied the evolution of bifurcation points 
and chaotic thresholds due to changes in h and /^bT, assuming always h ^ 0, ksT ^ 0. We 
have obtained the bifurcation point for the last Hopf transition and the last chaotic threshold 
for different values of ksT and h. In all cases it is found that both points are shifted towards 
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lower values of Vq with increase of either fc^T or h, agreeing with previous results j21|, l23j . 
which deduce this behavior from the reduction of quantization and temperature effects to an 
effective noise. This shift seems to have quadratic scaling in both parameters, in difference 
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This is not surprising taking into 



with the non-integer exponent change given in 
account the strong truncation in the power series we have performed. 

The temperature effects can be seen in Fig. (7), where we plot the value of the potential 
parameter Vq for which the systems enters the chaotic regime, as function of the temperature 
ksT. The different curves correspond to fixed values of h. The solid line set of curves 
correspond to the beginning of the Hopf bifurcation where A = 0, and the dashed ones are 
for the beginning of fully developed chaos. Note that even when this curves are not as 
smooth as the first ones, the general tendency is maintained for all of them. It can be seen 
how the value of Vq decreases as fc^T is raised, that is, keeping everything else constant, the 
greater the temperature, the earlier chaos appears in the system. The cutoff frequency has 
some effect on the curves, but the general behavior is the same independently of the value 
oi Uc- The influence of the cutoff frequency is more clearly seen in Fig. (8), which shows the 
same two transitions, for fixed ksT and variable h, and in Fig. (9) where it is clear that an 
increase in the cutoff has a similar effect as augmenting the h or the temperature, i. e., it 
shifts the thresholds for chaos towards lower values of Vq. The scaling of the threshold and 
bifurcation values with fc^T is quadratic as mentioned before. This can be seen in Fig. (10). 
This is probably due to the truncation to second order the whole EVM. 

A very similar behavior is obtained when we change the value of h. In Fig. (11), the 
evolution of bifurcation and threshold points is shown for several fixed values of ksT. The 
tendency as before is to have an early onset of chaos for larger values of h. However, a 
peculiar phenomenon appears here: for small values of the cutoff parameter, this tendency 
is reversed, and the system gets stabilized when h begins to grow, see Fig. (12). This 
anomalous behavior does not show up when we change ksT. 



V. CONCLUSIONS AND FINAL COMMENTS 



The kicked harmonic oscillator with dissipation has been quantized by means of the 
Quantum Langevin Equation formalism, truncating the moment expansion up to second 
order. We obtained the evolution maps for the expectation values Q and P and the second 



11 



order fluctuations i^P)"^ and SQSP + SPSQ. The most important result we flnd 

is that chaotic behavior appears in a region of parameter space far below the classical 
threshold for chaos. This agrees with the result reported by Liu and Schieve [3]; who also 
found chaotic behavior in a system whose classical version is periodic. We also find that the 
change in dimensionality in the map changes the type of bifurcations in the dynamics, as 
shown by the appearance of Hopf bifurcations on the quantized system. It should be noticed 
that this is not an artifact of the EVM, since any representation of a quantum system via 
c— numbers (wavefunctions in configuration space, elements of the density matrix in any 
given base, etc.) is infinite dimensional. Therefore, the reduction of a quantum dissipative 
system to the corresponding classical system plus some noise ignores the possible effects of 
this change in dimensionality. 

Once we are in a finite ksT and h domain, a continuous displacement of chaotic thresholds 
and bifurcation points towards lower values of Vq is observed, as fc^T and h are increased. 
This agrees with most the results obtained up to now in this field, and goes to confirm the 
rule that states that quantization of Hamiltonian systems eliminates chaos, hut dissipative 
systems become more chaotic when they are quantized. This is due to the loss of coherence 
brought in by the dissipative environment. By now, there are some experimental confirma- 
tions of this rule Here we should comment that irregularities observed in the curves 
for the threshold of chaos in Figs. (7) and (11) are due to the uncertainties incurred in the 
numerical evaluation of the largest Lyapunov exponent for the system. 

We are aware, of course, that a second order expansion is not sufficient to reproduce the 
behavior of the infinite dimensional EVM. Hoverer, we doubt that the inclusion or higher 
order corrections will restore periodic behavior to these new chaotic regions, although it is 
of course expected that some extra phenomenology will appear with the inclusion of these 
terms. 

We have introduced a cutoff frequency ujc in order to perform some of the integrals. These 
cutoff frequency is expected to come from the dimensions of the physical system in a real 
experiment. The existence of a cutoff implies a breakdown on the Markovian condition, 
but, as long as Uc ^ assuming that the system is Markovian is not a very strong 
approximation. In our simulations, we find that this is not a very important parameter 
except when the obviously inconsistent case when it takes very low values. 

Acknowledgments: M. R. wants to thank Cinvestav-IPN and the Fondo Yucatan for 
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Appendix A: Expectation Values Mapping 

If we are given the collection of moments 



= Tr(pp"a;"^) (1) 

= J dx{x\pp''x"'\x), (2) 

p being the density matrix of the system, we can construct the function 

Cia,P) = E (3) 

n,tlo-^ m!n! 

= j dx {x\p exp (iap) exp {iPx)\x). (5) 

This can also be written as 

Cia,P) = /da;cip(a;|pexp(.ap)|p)(p|exp(./3.)|a;), (6) 

= J dx dp (x|p|p)e'"^'+^^^-^*'^/^. (7) 

If we Fourier-transform the above equation we get 

j dad/3e-'''P'e-'^'''C{a,p) = j dx dp {x\p\p)e-'P''^H{p - p')S{x - x'), (8) 

(a;'|p|p')e-'^'^'/^ (9) 



47r2 



Solving for {x'\p\p') and changing the notation leads to 



{x\p\p) = ^^e^f^/'^y dad(3e-'''Pe-'^''C{a,(3). (10) 

If we want to get p in a position basis, we can take the left hand side of the above equation 
and do the following 

/ dx' {x\p\x'){x'\p) = [ dx' (x|p|x')e^P^'/^ (11) 

J V 2Tih J 

multiply both sides by e"'^'^^" 1^ and integrate over p 

= / cia;'(ip(x|p|x')e*f(^'-^")/^ = ^Ph^^,x\p\x") , (12) 
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and obtain the desired result 

{x\p\x") = dpe'P^''-''"^'^ j dadl3e-'^P-'^''C{a,l3). (13) 

This can also be put in the form 

{x\p\x') = ^jd(5C (^,/?) e-^^ (14) 

Appendix B: Expectation values of the inhomogeneous terms 

In order to calculate the 5 maps for Q, P and its three fluctuations, we will need the 
following terms 

{f-f)^f-f, (15) 

(hn-hri) = K-K, (16) 
if hr^ + hr^ f) = f-hn + hn-f. (17) 

We make the following definitions 

Is{t,uo,-f) = — dt'e-'^^*-^'^ sin ujo{t - t')F{t'), (18) 
Ic{t,u;on) = — t dt'e-^^'-''^ coscuoit -t')F{t'). (19) 

Uq J to 

The quantities we need are now expressed as 

f.f^^I^.I^, (20) 

UJq 

hn-hn= \ [(7 + ^"(Qn))'/. " + U^lh ' h ' Ml + V"{Qn)){Is ' h + h ' 



f ■ K + K ■ / = A [-2(7 + V"iQ^))Is ■ Is + -Ic + Ic- Q] , (22) 
and the integrals are explicitly given by 

Gssii^^UJon^T): (23) 

Gcc{(^,U;on:'T): (24) 

Gsc+cs{i^:U;on:r), (25) 



Is 


Is = 


2-fh 


POO 

/ dujuj coth 






TT 


Jo 


Ic 


Ic = 


2-fh 


POO 

/ duu coth 


TT 


Jo 


Ic + Ic 


Is = 


2^h 


poo 

/ dwuj coth 




TT 


Jo 



2kBT^ 

huj ^ 



where 



Gss(a;, a;o, 7, r) = / dx f dye '^^^^^^ cosujix — y)s\n.ujQXSva.ujQy, (26) 
Jo JO 

Gcc{u},LJo,'f,T) = dx dy e~"'^^~^^^ COS Lj{x — y) COS uqx COS LVoy, (27) 

^0 ^0 

Gsc+csii^^i^on^r) ^ dx dye"^^^"'"^^cosa;(x-y)sina;o(x + y). (28) 

Jo Jo 
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The integrals ()26ll28p are evaluated exactly, and give expressions that go like cu"^ for 
large u). With them, and the given values of ojq, 7, r, ksT and fi, integrals ()23H25p are done 
numerically, including a cutoff factor e'^^^" to control the logarithmic divergence. A Taylor 
expansion of coth is needed for c^o ^ 0. 
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FIGURE CAPTIONS 

Figure 1. Bifurcation diagram for the classical system for 1 < Vo < 9. Here we have set 
7 = 0.03, T = 10, and use 100 random initial conditions for each value of Vq. 
Figure 2. Attractor for the classical system using the same parameters as in Fig. (1), and a 
value of Vq = 8.0. 

Figure 3. Bifurcation diagram for the quantized system. We use the same values of 7 and r 
as in the classical system with a cutoff frequency cuc = 25. Here 100 random initial conditions 
are used, and h — ksT — 0.0002. Only the central part of the range in Q is shown for clarity. 
Figure 4. Magnification of the quantum bifurcation diagram in the second chaotic region, 
showing the presence of the Hopf bifurcation. 

Figure 5. Comparison of the largest Lyapunov exponents as function of V^. The dotted 
line corresponds to the classical case, which undergoes a period doubling bifurcation. In the 
quantum case (solid line) this is replaced by a Hopf bifurcation followed by a chaotic strip. 
Figure 6. Attractor for the quantum system using the same parameters as in Fig. (3), and a 
value of V^) of 8. Only the central part of the ranges in both P and Q are shown for clarity. 
Figure 7. Plot showing the shift of the chaotic thresholds (dotted hue) and the Hopf bifur- 
cation (solid line), with respect to ksT, for two values of the cutoff frequency. Each curve 
within any of the four bundles corresponds to a different value of h. Starting from the upper 
curve the values of h are 0.0005, 0.002, 0.004, 0.006, 0.008 and 0.01. 

Figure 8. Effect of the cutoff frequency on the chaotic threshold near Vq — 3.92. Values for 
ksT and h are fixed at 0.0001 and 0.0005, respectively. 

Figure 9. Beginning of the Hopf transition as a function of ksT for h — 0.006. with different 
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values for the cutoff frecuency, Uc — 1, 5, 25 and 125, witfi Uc augmenting downwards. 
Figure 10. Logarithmic plot of the beginning of the Hopf transition as function of ksT. The 
thin line has slope 2. 

Figure 11. Plot showing the shift of the chaotic thresholds (dotted line) and the Hopf 
bifurcation (solid line), with respect to h, for two values of the cutoff frequency. Each curve 
within any of the four bundles corresponds to a different value of ksT. Starting from the 
upper curve the values of ksT are 0.0005, 0.002, 0.004, 0.006, 0.008 and 0.01. 
Figure 12. Anomalous behavior of a chaotic threshold (upper bundle) and Hopf transition 
(lower bundle) for a low value of the cutoff frequency, as function of h. Starting from the 
upper curve the values of /c^T are 0.0005, 0.002, 0.004, 0.006, 0.008 and 0.01. 



17 



This figure "FIG01.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG02.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG03.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG04.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG05.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG06.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG07.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG08.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG09.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG10.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG11.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



This figure "FIG12.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/quant-ph/040400 1 v 1 



